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ABSTRACT 

The whispering gallery modes (WGMs) of optical resonators have prompted intensive research efforts due to 
their usefulness in the field of biological sensing, and their employment in nonlinear optics. While much infor¬ 
mation is available in the literature on numerical modeling of WGMs in microspheres, it remains a challenging 
task to be able to predict the emitted spectra of spherical microresonators. Here, we establish a customizable 
Finite-Difference Time-Domain (FDTD)-based approach to investigate the electromagnetic WGM spectrum of 
microspheres. The simulations are carried out in the vicinity of a dipole source rather than a typical plane-wave 
beam excitation, thus providing an effective analogue of the fluorescent dye or nanoparticle coatings used in 
experiment. The analysis of a single dipole source at different positions on the surface or inside a microsphere, 
serves to assess the relative efficiency of nearby radiating TE and TM modes, characterizing the profile of the 
spectrum. By varying the number, positions and alignments of the dipole sources, different excitation scenarios 
can be compared to analytic models, and to experimental results. The energy flux is collected via a nearby 
disk-shaped region. The resultant spectral profile shows a dependence on the configuration of the dipole sources. 
The power outcoupling can then be optimized for specific modes and wavelength regions. The development of 
such a computational tool can aid the preparation of optical sensors prior to fabrication, by preselecting desired 
the optical properties of the resonator. 
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1. INTRODUCTION 

Whispering gallery modes (WGMs) are produced by the propagation of electromagnetic waves along the interface 
of the surface of a resonator and its environment. As a result, an evanescent field is produced, which extends 
beyond the surface of the resonator. WGMs represent an optical phenomenon which is applicable to sensing 
applications, due to the sensitivity of the evanescent held to the presence of virions,^ and macromolecules such 
as proteinswhich break the symmetry of the held.l^ WGMs have been used for the development of label-free 
detection technologiespE optical frequency combs,!^ nonlinear opticJ^^and quantum electrodynamics (QED).I^ 
This work primarily focuses on the establishment of a computational modeling tool for investigating properties 
of WGM spectra in microspheres and microshells. 

Optical microresonators have been the subject of much recent investigation since the discovery of extremely 
high quality(Q)-factor generation (> 10®), which is possible by exciting the WGMs of a resonator, using a 
prism, or a tapered hber coupled to the resonator An alternative approach to this is to excite the WGMs 
indirectly, using fluorescent media, ^ surface plasmon effect^SEH or nanoparticle coating/doping using quantum 
dots.ll^ WGMs of the resonator also exhibit high sensitivity to their external environment. This means that not 
only are WGMs potentially very distinguishable from one another in the power spectrum, but that changes due 
to the presence of a new electromagnetic source can also be identihed easily. 
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The manufacture of optical resonators suited to a particular purpose using trial-and-error methods is expen¬ 
sive. Efforts to simulate and characterize the optical properties of WGMs in resonators using computational 
methods have so far focused on specific and ideal scenarios, such as excitation from a plane-wave beam.l^SHUl 
This opens the way for an investigation into realistic mode excitation scenarios, and analyzing differences in 
the spectral profiles, with a view to making a robust prediction tool. Here, we investigate coupling strategies 
that involve the excitation of WGMs of microresonators, using dipole sources placed on the surface. The Finite- 
Difference Time-Domain (FDTD) computational tool developed in this work will address this coupling scenario 
in particular. The tool is also easily extended to include fiber-coupled methods of mode excitation. 


2. FDTD MODEL 

The FDTD techniqu^^ splits a volume of space into a three-dimensional grid, and calculates the solution of 
Maxwell’s Equations for the electromagnetic fields at each point. This is repeated for a finite number of time in¬ 
crements, and performed sequentially as the system evolves in time. FDTD constitutes a comprehensive solution 
of an optical system, with all radiation modes present in the simulation. The only artifacts in the calculation 
involve discretization and finite-volume effects, and assumptions made regarding ideal material properties and 
resonator configurations. The FDTD method is also suitable for accessing transient or unexpected electromag¬ 
netic properties, as each time step is evaluated individually. However, the method is computationally expensive 
when achieving a result that is sufficiently converged.!^ 

The power spectrum represents a useful quantity for representing the positions and the Q-factors of the excited 
modes.l^ To obtain this, a flux-collection region may be introduced into the simulation, placed a certain distance 
from the microresonator. The total output power (P) in terms of wavelength (A) is obtained by integrating the 
projection of the Poynting vector (S = E x H*) onto a flux region of area A, with a normal unit vector n 

P(A)= [sndA. (1) 


The profile of the power spectrum, plotted as a function of A, will show sharp peaks that correspond to the 
positions of TE and TM WGMs. Ghanges in this energy output exactly map out the WGMs, and thus offer an 
excellent diagnostic tool for the quality of the resonator scenario being examined. 

An important aspect in developing the FDTD tool is a comparison with typical analytic models used in 
the fieldpHEZl which can provide estimates on the systematic uncertainty of the WGM positions. To check the 
mode positions predicted by FDTD, an analytic model is used, in which Maxwell’s Equations are solved in a 
dielectric medium with a boundary.!^ In this model, the resonance conditions for the TE and TM modes may 
be expressed in terms of the spherical Bessel functions of the first kind, j„i, and the spherical Hankel function, 
hm\ for azimuthal mode number m. These conditions are 

ni (m -K l)jm{z^) - z^jm+ijz^) ^ nz 
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respectively, for a refractive index of ni inside a sphere of radius R, and nz outside the sphere. The size 
parameters are defined as z ^2 = ^ ^^ 1 , 2 , for a wave number k = w/c, where c is the speed of light in a vacuum. 
Accurate matching of the mode structure of the spectrum allows one to identify the dominant contributions 
to the features of the spectrum. This identification facilitates future experimental studies in the reduction or 
enhancement of these features. 
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3. SIMULATION RESULTS 

The FDTD method will be used to consider a variety of flux collection scenarios, source distribution config¬ 
urations, and resonator properties. The use of dipole sources in the vicinity of a microresonator allows one 
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Figure 1. A circular flux collection region, offset from a 6 fim diameter microsphere, is placed so that Its normal is aligned 
perpendicular to the surface of the sphere. In this example illustration, the optical modes are excited from a tangentially 
oriented electric dipole source. 


to investigate the dependence of the optical modes on the size, shape and refractive index contrast associated 
with the resonator and its environment. Efforts to characterize the optical properties of WGMs in resonators 
using modeling techniques include both excitation from a plane-wave beam,t2®^^ and excitation from an electric 
dipole source.l^^lt^ The dipole sources can also serve as an effective analogue for fluorescent dyes or embedded 
nanoparticles that excite WGMs in microspheres 

3.1 Microspheres 

An FDTD simulation of a polystyrene (n = 1.59) microsphere, with a diameter of 6 /rm, is carried out, as 
illustrated in Figure [l] To generate an electromagnetic current to excite the WGMs, one or several electric 
dipole sources may be placed in the vicinity of the sphere. Here, an electric dipole source is used, which emits 
a Gaussian pulse with a central wavelength of 600 nm, and a width of 5 fs. In simulating the orientation of 
the dipole with respect to the surface of the sphere, both the tangential and radial cases are considered. Note 
that the pulse width is significantly narrower than the decay transition rate expected from a typical fluorescent 
source, such as Rhodamine dye, which is approximately 1-3 nsESIHl A circular flux collection region with a 
diameter of Influx = 2.58 fim is placed a distance of Tflux = 240 nm from the surface of the sphere, with its 
normal perpendicular to the surface. The power spectrum is then collected for wavelengths in the range 500-750 
nm. The spectrum is normalized to Pq, the dipole emission rate in an infinite bulk medium of the same refractive 
index as the surrounding medium. 

A comparison of the analytic model and the FDTD simulated spectrum for a microsphere is shown in Figurej^ 
for tangentially and radially oriented dipole sources in both air and water media. The WGM positions from the 
analytic model are shown as vertical bands, with TE^.n modes in green and TM^.n modes in red, for azimuthal 
and radial mode numbers m and n, respectively. The width of each band indicates the estimate of the shift 
in the WGM positions due to the uncertainty in the sphere diameter, 6 /im ±Aa;/2. In this case, the spatial 
resolution is held fixed at the value Aa: = 22 nm. For the temporal resolution, the spectral density yields an 
uncertainty of 0.31 nm. The two uncertainties are added in quadrature. 

3.2 Microshells 

The FDTD method can be modified to consider a microshell resonator, where a hollow sphere is placed in a 
medium, and excited from a number of dipole sources. The shell may be of variable thickness, and in this 
investigation, is made from a higher refractive index material than the external medium. The interior of the 
shell is assigned the same refractive index as the surrounding medium. To approximate the effect of exciting the 
WGMs using a fluorescent coating, an evenly-spaced distribution of electric dipole sources is simulated on the 
outer surface of the microshell. The orientation of the dipoles is randomly generated, and the number of dipoles 
placed on the surface (in this case, 75) is made large enough to provide an effective analogue of a medium of 
fluorescent emitters. Figure shows the variation of the power spectrum due to changes in the shell thickness. 












A [jum] A [/im] 

Figure 2. FDTD simulation of the normalized power spectrum of a polystyrene microsphere, 6 /rm in diameter, with a 
surrounding medium of air. Whispering gallery modes are excited from a tangential {top left panel) or radial {top right 
panel) electric dipole source with a central wavelength of 0.6 fim. The results for a surrounding medium of water are 
also shown for a tangential source {bottom left panel) and a radial source {bottom right panel). Vert ical lines indicate 
predictions of the TEm,n {green) and TMm,n {red) modes derived from a typical analytic modelp2I2Il for azimuthal and 
radial mode numbers m and n, respectively. The width of the bands indicates the systematic uncertainty in the positions 
due to the finite grid size of FDTD. 




Figure 3. Power spectrum simulations of fluorescent microshells with a shell refractive index of n = 1.45, and with a water 
medium {n = 1.33) both inside and outside the shell. The modes of the shells are stimulated using an evenly-spaced 
distribution of 75 electric dipole sources of random orientation, placed on the outer surface of the shell. Left panel: the 
outer-shell diameter is 10 fim with shell thicknesses of 60, 75 and 90 nm. Right panel: the outer-shell diameter is 20 ^m 
with shell thicknesses of 120, 150 and 180 nm. A circular collection region (diameter: 2.58 pm) is placed 240 nm from 
the surface of the shell. 







































































































































































































In the left panel, an outer shell diameter of 10 /rm is used, with shell thicknesses of 60, 75 and 90 nm considered. 
In the right panel, an outer shell diameter of 20 /rm is used, considering shell thicknesses of 120, 150 and 180 
nm. For a low refractive index contrast, a small sphere diameter and a small shell thickness, the WGMs are 
not easily resolvable. However, as the shell thickness increases, an improvement in definition is observed in the 
spectrum structure. For spectra obtained from a larger shell diameter (shown in the right panel of Figure a 
peak structure begins to emerge as the shell thickness is increased. Note that the presence of the distribution of 
dipole sources serves to broaden the WGMs peaks by exciting a range of closely-spaced higher order modes. 

4. SUMMARY 

A customizable FDTD tool is established, and investigated in the context of whispering gallery mode generation 
in microresonators. The method can be tailored to more realistic scenarios, where the assumptions built into 
typical analytic models make direct comparison with experiment difficult. An important novel feature of FDTD 
is the ability to investigate the choice of source used for mode excitation, as-yet unexplored in the literature. 
Dipole sources may be placed at a variety of positions and orientations on the surface of the sphere, which can 
serve as an analogue for nanoparticle coatings. 

The generation of whispering gallery modes in novel resonator structures represents an important applica¬ 
tion of the customizable FDTD approach. The FDTD tool is extended to explore configurations of microshell 
resonators, represented by a hollow sphere with a water medium placed both inside and outside, where the shell 
is made from a higher refractive index material. A distribution of randomly-oriented dipole sources is placed on 
the surface, analogous to a fluorescent dye coating. 

The FDTD tool presented here represents the first step in establishing a realistic optical configuration sim¬ 
ulator, useful for facilitating a cost-effective approach to designing tailored optical resonators. Tuning the 
characteristics of the resonator for optimal power coupling and identifying the most prominent spectral features 
will help to reduce fabrication costs and aid the development of the next generation of biosensing tools. 
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